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FOREWORD 

The researches described here were supported mainly by the National 
Aeronautics and Space Administration (under Grant NCR 05-007-416) and by 
the UCLA School of Engineering and Applied Science. These studies were 
part of a continuing program of researches In gas-surface Interactions. 
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ABSTRACT 


Energy transfer In collisions of satellite-speed (7000 m/sec) helium 
atoms with a cleaned 6061-T6 satellite-type aluminum surface was investi- 
gated using the molecular-beam technique. The amount of energy transferred 
w.is determined from the measured energy of the molecular-beam and the 
measured spatial and energy distributions of the reflected atoms. 

Spatial distributions of helium atoms scattered from a 6061-T6 aluminum 
surface were measured again in this study, and show features similar to 
those presented in report UCLA-ENG-7546 [1]. The scattering pattern exhibits 
a prominent backscattering, probably due to the gross surface rougnesss and/ 
or the relative lattice softness of the aluminum surface. 

Energy distributions of reflected helium atoms from the same surface 
were measured for six different Incidence angles. For each incidence angle, 
distributions were measured at approximately sixty scattering positions. 

At a given scattering position, the energy spectra of the reflected helium 
atoms and the background gas were obtained using the retarding-f ield energy 
analyzer. The mean ref lected-beam energy and the differential energy accom- 
modation coefficient ((A.C. )_(0. ,6 ,4>)) were then extracted from these 

b 1 r 

spectra using a least-square fitting program. The measured (A.C.) (6,,0 ,i}))s 

b 1 r 

show some fluctuations and a weak dependence on scattering angle, i.e., the 
accommodation decreases slowly as the scattering direction shifts toward 
the surface tangent. 

The overall energy accommodation coefficient for a beam with a given 
Incidence angle was then evaluated using the measured spatial density dis- 
tributions and the mean ref lected-beam energy distributions. Results show 
that the mean accommodation coefficient varies between SOX and 65%, depen- 
dent on the incidence angle. 
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INTRODUCTION 

Banlc knowledm.' concerning energy and nio.. ntutn transfer between earth 
satellites iiid upper-atmospheric gases is essential for understanding the 
drag experienced by earth satellites (therefore for estlnuitlng the life- 
time of an earth satellite and/or extracting the mean upper-atmosphere 
dc'islty from satellite drag data). For example. In predicting the aero- 
dyn.j.Tilc drag of a satellite, one uses frequently a model In which the 
thermal accommodation between the ambient gas and the satellite surface 
is complete and the scattering distribution of reflected molecules follows 
the cosine law. However, possible deviations from this model might yield 
greatly different results. 

These energy and momentum transfers can be Investigated experimentally 
In the laboratory using an ultra-high vacuum system and the molecular-beam 
technique. The desired Information can be extracted from the change In the 
beam properties during the surface collision If the states of both the inci- 
dent and the scattered beam (spatial distribution and speed distribution) 
can be determined. Spatial distributions of satellite-speed helium beams 
scattered from satellite surfaces wen* obtained previously and summarized 
In report UCLA-F.NG-7546. This report presents measured energy distributions 
of helium atoms reflected from 6061-T6 aluminum surfaces. 

In Chapter 11, the experimental apparatus and procedures are described 
briefly. Emphasis Is given on the design and the operating procedure for 
the retard Ing-f leld energy analyzer. Experimental results are given and dis- 
cussed In Chapter III. A least-square curve-fitting computer program Is 
given In an Appendix. 
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CHAPTER II 


KXPERIMENTA'. APPARATUS AND PROCEDURES 

The present experimental study was carried out In the UCLA Molecular- 
Beam Laboratory using the molecular-beam system shown schematically in 
Figure II-l. Since it has been described In detail elsewhere (1,2], o,ily 
a brief description will be given here. 

The satellite-speed (7000 m/sec) helium beams were generated using an 
arc-heated supersonic beam source developed by Young (3]. The Incident 
beam was. collimated by an orifice of 0.10-lnch diameter placed between 
the colllmatlon chamber and the detection chamber. The beam was character- 
ized by a multi-disk velocity selector located in the colllmatlon chamber. 

A n>w detection system was constructed during the course of this study 
for facilitating measurements of the complete three-dimensional density and 
mean-energy distributions of satellite-speed helium atoms reflected from 
satellite surfaces. Cf. Figure 11-2. This new system includes (1) a target 
positioning mechanism, (2) a detector rotating mechanism and (3) a mass 
spectrometer and/or a retarding-f leld energy analyzer. Descriptions of the 
first two mechanisms were given in the first report of this study (cf. ref. 1 
The design and the operating procedure for the retardlng-f leld energy ana- 
lyzer will be given here. 

The retarding-f ield energy analyzer is shown in Figures II-3 and II-4. 

An electron-impact ionizer, mounted 0.5)-inch from the target s irface on the 
entrance plate of the analyzer, was used to ionize a fraction of the beam 
species (also of the residual background). The retaiding-f d section of 
the analyzer is made of seven thin stainless-steel washer-shaped discs 
placed in a stainless-steel can. T' *nlet plate is followed by three 
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Figure IM. Schemat. Diagram of the Molecular Beam Syitem. 
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Figure 11-4. Schematic Diagram of the Retarding-Fie<d Energy Analyzer. 



tocusln^; plates, a retarding plate and two more focusing plates. The 
potentials of all plates except that of the Inlet plate are floated rela- 
tive to the retarding potential so that Ions of different energies will 
experience the same focusing effects within the analyzer. Thus undesir- 
able effects of the energy-dependence of the transmission efficiency 
are minimized. Typical plate potentials also are given In Figure II-4. 

The ions that have passed through the retarding-f ield region were 
filtered by a 2-inch quadrupole mass filter to eliminate the noise from 
the Ionized background gases. The filtered ions were then detected by a 
p Ise-countlng particle detector. The energy spectrum of the reflected 
atoms at a given scattering, position was obtained by measuring the 
reflected-beara density as s. function of the retarding potential. The 
measured spectrum was proct;ssed by a NS513 signal averager and recorded 
on IBM cards. A block diagram of the electronic system is shown in Figure 
11-5. 

Although the elec tron-lmpact ionization does not change the kinetic 
energy or a helium atom (since the translational energy transfer between 
the lonizlr-; electron and the atom is negligible due to the large ratio of 
their masses), it was found that space-charge effects of ionizing electrons 
In the ionization region and/or surface-charge effects on the anode cage 
did Introduce a systematic shift of the entire energy spectrum toward lower 
energies (l.e., the positive ions were produced in a region of negative 
potential with respect to ground). To reduce this shift, a small emission 
current 50 uA) was used in the ionizer. Also, a positive potential (8 
volts relative to ground) was applied to the anode cage in order to counter 
shift the energy spectrum toward higher energies. Then, since the potential 
of the ionization region was no longer at ground level, it was necessary to 
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Figure 11-5 Block Diagram of the Electronic System of the Retarding-Field Energy Analyzer. 















nscercain Independently a reference point in the energy spectrum. This 
reference point was provided by the thermal energy spectrum of the back- 
ground gas, which has a mean thermal energy of 0.05 eV (l.e., the mean 
thermal energy at 296*K). 

Since the background gas of the beam species also contributed to the 
measured spectrum, it was necessary to subtract this contribution in order 
to obtain the ref lected-beara energy distributions. This subtraction was 
facilitated by measuring two spectra (one for the reflected beam plus back- 
ground and one for the background alone) under the same operating conditions. 
Both spectra were then processed using a computer program; the reflected- 
beam energy spectrum was obtained by subtracting the background spectrum 
from the overall beam-energy spectrum. Both the background spectrum and 
the reflected-beam spectrum were least-square fitted using .. high-order 
Chebyshev polynomial function. The differential energy distributions f(E) 
were obtained by simple differentiation of the fitted functions. The mean 
reflected-beam energy at a given scattering position was evaluated from 

E (0, ,6 ,4>) = E - E , + 0.05 (eV) (II-l) 

r 1 r r ref 


where 


E(e^.6j.,<|)) = y*f(E)-E-dEyy*f(E)-dE 


(n-2) 


and 0.05 eV is the thermal energy of the background gas at 296®K. The dif- 
ferential energy accommodation coefficient at a given scattering position 
was obtained using 


[A.C.]^(e^,0^,r{)) = 


Ei-Er(ei.0r.4)) 

E^ 


(II-3) 
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where is the incident-beam energy. The computer program and Its input 
parameters for handling the described data reduction are given in the 
Appendix. The overall energy accommodation for a given incidence angle 
was then evaluated by 

[A.C.]g(6^) = £ £n^(e^,e^,4>)-[A.C.]g(0^,6^,<(.) (ir-4) 

0 4 ) 

r 

where n^ (0^ ,6^,$) is the normalized spatial density-distribution function 
of reflected helium atoms. 

As Indicated, spatial distributions of satellite-speed helium beams 
scattered from four different satellite surfaces were obtained in the first 
phase of this study. Experimental procedures and results are included in 
report UCLA-ENG-7446 [1]. However, spatial distributions at some angles in 
the backscattering region were not measured at that time due to the con- 
straint on the detector path as Indicated in Figure II-6-(a). This problem 
was solved later by rotating the surface counterclockwise beyond the normal 
incidence angle while retaining the previous detector path as shown in 
Figure II-6-(b). Spatial and energy distributions were measured in the 
present study using these complementary configurations. 
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DETECTOR PATH 

(IN THE PLANE OF INCIDENCE) 


CONFIGURATION (a) 



CONFIGURATION (b) 


Figure 11-6. Complementary Beam-Surface Configurations of the Scattering System. 
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CHAPTER III 


RESULTS AND DISCUSSIONS 

Spatial distributions of a satellite-speed (7000 m/sec) helium beam 
scattered from a cleaned 6061-T6 aluminum satellite surface for six dif- 
ferent incidence angles (0*, 15®, 30®, 45®, 60® and 75® from the surface 
normal) are shown In Figures IIl-l to I1I-6. The center of the polar dia- 
gram corresponds to the point of Impingement. The Incident beam Impinges 
on tlie test surface (which coincides with the surface of the page) from 
the bottom of the diagram with the given Incidence angle measured from the 
surface normal. The upper (6^ > 0) and lower (0^ < 0) halves of the dia- 
gram represent the forward-scattering and backward-scattering regions 
respectively. The dashed lines at constant value of 6^ indicate detector 
paths (l.e., from ■ 0® to (f> - 90®). 0 denotes the out-of-plane scatter- 

ing angle and 4* “ 0® represents the plane of Incidence. These results show 
diffusive scattering patterns and exhibit trends similar to those previously 
reported [1]. As Indicated before, the most Interesting feature ou these 
scattering patterns Is the prominent backscatter Ing of the Incident helium 
atoms (l.e., a large fraction of the Incident atoms are scattered back In 
the vicinity of the incident beam), particularly as the Incidence angle 
Increases toward the surface tangent (l.e., for large values of 6j^)» This 
large fraction of backscatter Ing could be due to the gross surface roughness 
and/or the relative lattice softness of the aluminum satellite surfaces. 
Smith [4] observed a large Increase In backscatterlng Intensity for Increas- 
ing surface roughness In his computer simulation of gas molecule reflections 
from rough surfaces. This backscatterlng could result in relatively high 
drag coefficients for such satellite surfaces. The spatial-dlstrlbutlon 
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Figure IIM. Polar Plot of Scattered -Beam Daniity Distribution for 7000 m/sec 
Helium Beam Scattered from Cleaned 6061 -T6 Aluminum Plate at 
0° Incidence Angle. 
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Figure III-2. Polar Plot of Scattarad Baam Dantity Diitribution for 7000 m/«ac 
Halium Beam Scattarad from Claanad 6061 -T6 Aluminum Plata at 
15** Incidanca Angle. 
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Figure MI-3. Polar Plo. of Scattarad-Beam Density Distribution for 7000 m/sec 
Halium 3aam Scattered from Cleaned 6061 -T6 Aluminum Plata at 
30° Incidence Angle. 
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Figure 111-4. Polar Plot of Scattered -Bean Deniity Diitribution for 7000 m/$ec 
Helium Beam Scattered from Cleaned 6061-T6 Aluminum Plate at 
45** Incidence Angle. 
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Figure III 5. Polar Plot of Scatterad-Baam Daniity Oittribution for 7000 m/tac 
Halium Baam Scattarad from Claanad 6061 -T6 Aluminum Plata at 
60° Incidanca Angle 
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Figure 1 11-6. Polar Plot of Scattered-Beam Density Distribution for 7000 m/sec 
Helium Beam Scattered from Cleaned 6061 -T6 Aluminum Plate at 
Tj° Incidence Angle. 
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meaKuremPnts shown here and the energy-distribution measurements to be pre- 
sented next provide the data required for estimating the overall energy 
accommodation for this beam-surface combinatiin. 

Measurements of energy distributions of satellite-speed helium atoms 
scattered from a cleaned 6061-T6 aluminum satellite surface were made for 
six different Incidence angles (6^ ■ 0“ , 15®, 30®, 45®, 60® and 75* from 
the surface normal). For each incidence angle, distributions were measured 
at approximately sixty scattering positions. These scattering positions 
included eleven In-plane scattering angles (6^ ■ ±75®, ±60®, ±45®, ±30®, 

±15® and 0®) and six out-of-plane scattering angles (4- » 0®, 15®, 30®, 45®, 
60® and 75®). Typical energy spectra obtained at a given scattering posi- 
tion are shown in Figures II1-7 and 1II-8. Curve A of Figure III-7 repre- 
sents the energy spectrum of a reflected helium beam superimposed on an 
energy spectrum of the background helium gas. Curve B of Figure III-7 
represents the (thermal) energy of the background helium gas (mostly due 
to beam load). The reflected :eam energy spectrum is then the difference 
of these two spectra (i.e., A-B). Figure III-8 shows the normalized energy 
spectra of the thermal background and the reflected helium atoms (Curves 3 
and 1), their least-square fitted curves (Curves 4 and 2) and the correspond- 
ing dif f'^rential energy distributions (Curves G and A) obtained using the 
computer program shown in Appendix A. The differential energy accommodation 
coefficient was obtained using Equations (II-l) and (II-3). Results for 
(A.C.)„(6.,6 ,4>) obtained at all possible scattering angles are given in 

£< X IT 

Tables III-l to III-6. Measurements v. re not possible wit^'in a solid angle 
around the incident beam (due to interference between th» detector and the 
incident beam at these scattering positions) and for some glancing scattering 
angles (due to weak signal-to-noise ratios). These tables also Include 
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RFTAPPING POTENTIAL (VOLTS) 

EIonkF III-7 ENFRGY SPFCTRA OF T'lE REFLECTED HELIUM ATOMS 
SUPERIMPOSED ON THE BACKGROUND (CURVE A) AND 
THF HERMAL BACKGROUND HELIUM GAS (CURVE B) 
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TabI* MI-1. The Diffarantial Enargy Accommodation Coafficiantt and tha Normalized Spatial 

Oantity Distribution for 7000 m/sac Helium Beam Scattered From Qeaned 6061 -T6 
Aluminum Plate at 0° Incidence Angle 



(a) The Differential Accommodation Coefficient (%) 

(b) The Normalized Spatial Denrity (%) 

Standard Deviaiton (%) 
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Tabit III-2. Tha Differential Energy Accomnwdation Coefficients and the Normal<«.:;* Spatial 

Daniity Distribution for 7000 m/sac Helium Beam Scattarad from oiaanad 6061 -T6 
Aluminum Plata at 15° Incidence Angi« 



NOTE: 

(a) The differential Accommodation Coefficient (%) 

(b) The Normalized Spatial Density (%) 

(c) Standard Deviation (%) 
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TabI* MI-3. Tbt Diffarantial Enargy Accommodation Coafficiants and tha Normalizad Spatial 

Dantity Diitribution for 7000 m/iac Halium Baam Scattarad from Qaanad 6061 -T6 
Aluminum Plata at 30° Incidanca Anf^a. 



NOTE: 


(a) The Differential Accommodation Coefficient (%) 

(b) The Normalized Spatial Density (%) 

(c) Standard Deviation (%) 
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TabU III-4. Tht Differential Energy Accommodation Coafficientt and the Normalized Spatial 

Deniity Distribution for 7000 m/tac Helium Beam Scattered from Qeaned 6061 -T6 
Aluminum Plata at 45° Incidence Angle 



NOTE: 

(a) The Differential Accommodation Coefficient (%) 

(b) The Normalized Spatial Density (%) 

(c) Standard Deviaition (%) 
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TabU lll-S. Tha Diffarantial Enargy Accomodation Coafficiano and tha Normalizad Spatial 

Damity Diitribution for 7000 m/tac Halium Baam Scattarad from Oaanad 6061 T6 
Aluminum Plata at 60° Incidanca Angle. 



NOTE: 

(a) The Differential Accommodation Coefficient (%) 

(b) The Normalized Spatial Density (%) 

(c) Standard Deviation (%) 
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Table III-6. The Differential Energy Accommodation Coefficient! and the Normalized Spatial 

Dantity Diitribution for 7000 m/tec Helium Beam Scattered from Qeaned 6061 -T6 
Aluminum Plate at 75°lncidenca Angle. 



NOTE: 


(a) The Differential Accommodation Coefficient (%) 

(b) The Normalized Spatial Density (%) 

(c) Standard Deviation (%) 
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standard deviations (o) of the reflected-beam energy-spectruo data from the 
least-square fitted curves and the normalized spatlal-dlstrlbutlon function 
of the reflected helium atoms obtained from the measured spatial distribu- 
tions shown In Figures III-l to III-6. The overall energy accomnodatlon 
coefficients at a given Incidence angle was then evaluated using Equations 
(11-4) and the data given in these tables. The results are shown in Figure 
III-9. 

The differential accommodations obtained show some fluctuations, due 
perhaps to the weak slgnal-to-nolse ratio which results from the relatively 
diffusive scattering from the satellite-type aluminum surface. The results 
also Indicate a weak dependence of accommodation on scattering angle, l.e., 
the (A.C. )g(6^,6^,0) decreases as the scattering direction shifts toward 
the surface tangent. 

The overall accommodation coefficient Is slightly higher for a glancing 
Incident beam than for a normal Incident beam. The value varies between 
50% and 65% for this beam-surface combination. 
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Figure III-9. 


Overall Energy Accomnrodetion Coefficient of a Satellite-Speed Helium Beam 
(1.02 eV) Scattered From a Cleaned 6061-T6 Aluminum Surface at a Function 
of the Incidence Angle. 
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APPENDIX 


COMPUTER PROGRAM FOR LEAST-SQUARE FITTING 


A. PROGRAM : 


C NAlN INMI.T; call; output 

c 

t I McNii I NC ( * I i ) t 

I t ( « I 2 » . Yew ( <^l 2 ) . Vaol t>l ^ ) t 

1 DYLS(tl2). 

1 CH2(2 ) .SIGWA2 (2 ) . 

t OF ( 2£ 1 t 1 0 1 ) .OV ( 2£) I 

CLWM,,\ /AAA/X(3l2)tVN(tl2)tYLS(bl2)tOY(512)«N.fc.M.SIuMA«KM 
CATA »JuANK./1m /.UOT/IM,/ 

Data e a/ jh i / , bml / i m2 / ,CFw/ iHd/ 

DATA Uii/lMj/,b * / I MA / ,Ctie/ I MC./ 

I CONTINOC 

WCAO (S.JOO) Nil T , XMCNtACA Y.xYLAt^t THLTA I .TMETAW.PmI 

IF (nsct.lt.i) go to lOCO 

•HITE (6.^00) KMuN tKC A V ,K YE A.3 , TMtT AI , T efcT/.H . PMI 
If- INSH»^«E0*2) GO TO E 
IF (NS>eT.uT,2) NiET = 2 

OFAD (S.FIO) r I . M, KC V I , V 1 , NC V2 .V2 .NI tXLL .Wr.XHL. 

1 INT I, INTF .ikPGtKSTEP.NPK ,NSTtP 

N sNF - N I ♦ I 

GV=(V2-V1 )/( NCV2-NCV I ♦ 1 I 
VI -VI aDv*! M -nCVI } 

VF = Vl*t)V*( NF-KCV 1 ) 
b CONT 1 NOE 

If ( KPG.LT. 1 I (jL to T6 

CO 10 f=1.251 ■ ' 

GV( I ) = 3 
CO IC J = I # lOl 
Cf ( I . J ) ANn 

10 CCNTInoE - -- - - 

DO 2C I * 1 . 25 I .*i 
CO 20 J= 1 . 1C I 
uF I I , J ) =l;ct 

3 0 CONTI NUF. ■ - . 

CO 30 I = I . 25 1 
DO 30 J=l . ICl .1C 
GF( I , j ) =UJT 

30 CCnTINuF ■ _ 

35 continue 

CO 40 1=1.32 

K = 1 e>*< I -1 I 

40 PEAD (5.230 C VhG(K«- J) ,j = l .161 
IF (NSET.LT. 2) GO TO 55 
CO SC 1=1.32 
1 5 •( I - I I 

READ (S.mO) (VFM(k ♦ J > . J = 1 . 16) 

DO bC J = 1 . 1 6 

•30 YiJM ( K ♦ J ) - YH4 ( K ♦ J ) - VO(j ( K ♦ J ) 
bS CONTINUE 

CO SCO KI=1.NSFT 
no 60 I =M . Nt 

F = I -M ♦ 1 
NC ( K I = I 
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nrr. 


i 


e 1 

t2 

AO 

70 


75 

80 


1 10 

1 20 
130 


1 00 
200 


201 


204 


202 


203 

300 


OU TCi (61 to2 1 .K 1 

V(K ) =VBG( I » 

CU TJ 60 
V( K I ZYHM( n 
CONTI NOC 
SYMAXzO .0 
SYMlNsO.O 
CO 70 1=1,50 
SYMlN = SY»>IN*t (►:-! » 
cyMAx=SYmAX4Y( I ) 

YMI N = S YMIN/50.0 
YMAX =0 Y MAX/50 *0 
DU 80 X = 1 • N 

YN(K)=(y(K)-YHINI/(YMAX-YMIN) 
IF (YN(K).CE.0*C I Ol, TC 75 
YN( K ) = 0* 0 

IF ( YN( K ) .Lb. I .0 ) OU TO 80 
VM K) = l .0 
CUNT INUE 


call LE A5T-SUtAFE CUUvE FITTINO 


CALL L SF Cm ( M , XLL. XHL ) 

GO TU ( 1 1 0, 1 20) ,K 1 
«R1 TE ( 6. )10) 

GO TO 130 
«»HITL (j.015) 

CUNT INUE 

IF (NPW.LT.I) go TO 200 
#<RTTF (6,020 ) 

CO 100 K=l,N,NSTEP 

•MITE (O.O30) NC ( K ) . X ( K ) , Y (K. ) • YN ( X ) , YL S ( K ) , OY( K ) 
CONTI NUE 

CONTINUE - — 

£M=0« 0 

sigma =0.0 

call E mean! I NT I , INTF ) 

EM2(KI ) =EM — — - 

call oevia 

b I GMA2 ( K ( ) = S IGMA 
51 GMAPsSIGMA *100.0 

W9ITE le.'OTO) cm.sigwap - 

IF (KPG.LT.1) GU to 400 
K I N= I 

CU 290 K=1.N,KSTEF 

KIN=KI NFl - 

GV( KIN ) = X( K ) 

JYN=100*AdS(YN(K) )♦! 

IF (JYN.LT.101) GO TC 291 

JYN= 10 1 • 

CONT 1 NUE 

IF ( YLS( K ) .GE.0.0 ) GC TO 2 94 
JYL 5= I 


GJ TU 292 
CONT I NUE 

JYLS=1 3 0*(A3S<YLS(t<) ) )♦! 

IF ( JYLS.LT.lOl ) GO TO 292 
JYL 5=1 01 

continue 

JDY = 10 0 *AHS( 0V( K ) ) ♦! 

IF (JUY.lt. 101) GO TC 293 
JDV = 10 1 
CCNT 1 NUE 

GO TfJ ( JC 0. 1 1 0 ) ,K I 
GF ( K 1 N , J YN ) =oG 
OF ( K IN, JYLS)=dGL 
GF ( K 1 N , JCY ) =DUG 


Of“&/-'QE IB 

QTJa j 


fiUAXJTY 
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non 


GO TO " ■ 

.3 10 Of ( K I N . JYM 

Gf < K IN , JVL S > =fiMt 
CF ( K IN • JCV ) »COM 

29“ CONTiNue '* 

2 JO continue 

IF ( NSE r «uT . 1 • ANO.K I ,LT .2 ) Gu TJ 400 
•PITE (6.<S40) 

GV( I IsGVC 2 )-( GV( 1»-GV{ ? ) ) 

CO 3e>0 ! *l « K I N 

4PITE (6*950) oV( I ). (OF ( I , J) *jal . 1 01 ) 

360 continue 

400 CONTINUE * 

500 CONTINUE 

CALCULATE THOL MfcAN-tJEAM ENEPGY 

ACCsO.O 

IF (NSET^LT.2) go to 600 
ETPUE»EM2( 2 )-EM2( 1 )4C«0» 

ACC = (EI-ETRUE )/EI*100 ~ 

600 CONTINUE 

(•RITE (6. J75 ) 

MRITE (6*960) ETRUE*ACC 

WRITE (6*975) “ 

GO TO 1 

dCO FORMAT ( 4 1 10 * 3F 1 0.3 ) 

810 FORMAT (F5.2* 15*4( 13*F7.2) *615) 

030 FORMAT (16F5.0) * ' " 

910 format (12X.23H** SET i: HACKCRJUNC *«/) 

930 FORMAT (// /2X * 7HDAT E : . I 2 * 1 H/ * I 2 * IH/ * I 2 * 5 X * 

1 2 6HANGLES( THET a I/TF-ETAR/PM 1 ) : *F6.2 * IM/ *F6. 2* 1M/*F6«2/ ) 

915 FORMAT (12X*17H*4 SET 21 BEAM •*/) “ * 

020 FORMAT (3X*7FCH. NO* . 1 2X . 3HV -R * 1 1 X * 5H 1 - S 1 G * 1 3X . 5H 1 -SOR * 

1 1 OX * 7ml SI -NUR *8X *9HC( LS 1 )/OV) 

9 30 FORMAT ( 5X * I 5 * 1 C X * F6 * 3 * 1 CX *F 6* 1 * 3 ( 10 X * F 7 * 4 ) I 

940 FORMAT ( // 8 X * 3 MV /R . 1 0 X * 7 H 1 : I ( UM ) * 3X* 9K2 TL S I ( BM I * 3X~* 

1 loHd :0( LSI ) /O V-( ) *3 X *7H3 : 1( UG ) * 3X *9h4 (LS 1 (PG ) * 

? 3X* 1 6MC (L5 I ) /OV- ( EG )// ) 

950 FCRMAT (2X *F 1 C*5 * 3X* 1 01 A1 ) 

J70 FORMAT (/5X , 1 2MVtAN-ENERGY : *F7*4 *3H EV , 1 0 X , 1 9 MST ANDARU DEVIATION:* 
1 2X*F5.2* IMX/) 

975 FORMAT! //lOX* 80 H****<<*4*4444***44 ** ** 4 * *» *4 ** *♦** *»-*^« (T* 

t ««*««•* ««*«**«•»«**«•«»«** «•««//) 

980 FORMAT (/lOX* 22MTRUE MEAN-BEAM EN ER GY : * 2 X *F 7 , 4 , 3X * 2 HE V * 1 0 X , 

1 26MACCOMMJOAT I CN CO EFF 1 C 1 EN T= • F 7 * 4 * 1 H «/ ) 

1000 CONTINUE 
STOP 
END 


SUBHOUTINE L SFCmv ( M* XLL * xml ) 

CCMMON*/AA A/X? 5 I 2 ) !y( 512 )*Yl 3 ( 512 )*OY( 512 ) *N *EM* SI GMA *KM 
COMMON /SIM/A( 100 ) *P( 10 ) _ 

Ml=MFl 

IM T I AL I ZAT ICN 
DO 2 I » 1 iMl 

?R(I) = 0*0 - 

I *2 = M 1 *M 1 
DO 4 1 = 1 *M 2 
It A( I ) = 0 *0 

XU = 2.0/(X(N)-X( I ) ) . - 

X 0 *(X(N) 4 X (1 ) )/(X(M-X(l )) 

CALCULATE' CREBYSmEV FUNCTION AT TX 
T( 1 )=t .0 

DO JO K = 1 *N _ ... . 

TX=XO*X|K 1 -XO i , 

T( 2 )=TX Cp iAf]!.} 


LaUIWI — - ^ 

T( I )=2*0*TX*T( I-l )-T { l-e ) 
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r n 


I 


03 2C I -1 . Ml 

P( I >=»( 1 )♦!( I >*V(K ) 

DO 20 jsl .Ml 
IJ=( 1-1 »•«!♦ j 
20 A( IJ) = A( IJ )♦!( r )*T( J ) 

30 continue: 

EPS=1 . OE-bO 
IER=0 

CALL SIMC(M1 ,1 .EPS. lEP) 

CALCULATL CEBbYSHEV MCLYNCMIAL AT TX USlNo rt(l) 
calculate OY FOPM CMEEYSEEV POcVNOMIAL 
0T( 1 ) =0. 0 
CT<2)=1.0 
DO 30 K = 1 . N 
YL?K = 0« 0 
nvK=o *0 
Tx*xo*x(*<:>-xo 

T( 2 ) = T X 
CO 35 I-J.Ml 

3£ T< 1 ) =2 . C»T X* T( I -1 1 -T ( I -2 ) 

DO 40 1=1. Ml 
40 YLSK= YLSK4T ( I ) ( I 1 
?0 YLS(K)=YLSK 
DO 60 1=3. Ml 

60 CT( I ) = 2.04(T( I-l )ATX*CT( I-l) 1-DT < 1-2 I 
CO 70 1=2. Ml 

TO CYK=OYK4P( I 1 AUT < I ) 

IF (OYK.LE.O.O) (iO TC 80 

DVK = 0. CCCOOOl ■ *• 

aO CV(K1=0YK 
DYM =0. 0 
CU 250 K=1.N 

IF ( X( K ) ,LT • XLL*OH .X< K ) ,GT«XHL 1 GO TO 250 
240 ’*= ( OY ( K ) . GT ,0 #0 1 GO TC 250 

Ir (OYM.GT.AaS(Oy(k) 1 ) GC TO 250 
CYM= AbS( OY ( K ) ) 

2£0 CONTINUE 
KRS=0 
K=0 

280 K=KAl 

KR = KM-Ka1 

IF (KR.LT.l) GO TO 200 
IF IKRS.GT*!) GO TO 265 
CY ( KH ) =OY ( KR 1 /OYM 
IF (DY ( KP) «L T.0«0 ) GC TO 280 
kRS = 2 - - - -- 

265 CY(KR) =-0,0C0001/OYM 
GO TO 280 
290 K=0 

KFS=0 ‘ - .. 

292 K=KF1 

KF = K M F K 

IF (KF.GT.N) GO TO 300 

IF (kFS.GT, 1 ) GO TO 2 94 ~ ' ’ “ 

CY ( KF ) =OY( KF l/OYM 

IF (OY ( KF 1 .LT.O *0 ) GC TO 292 

KFS = 2 

294 CV(KF»=-0,OC0C01/OYM 
GO TO 292 
300 continue 
RETURN 

END “ ■ ■ " 


“ ' JT. 
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1 


SUBtOUT I Nt Cfc V I A 

COM MOM /AAA/A(‘ili?» tVN(f'l2» .YLb(-il2)tCy(5»2 )tN.EM,SIOMA,KM 

SYrO.O 

SU = 0« 0 

DO 100 K = 1 tN 
SV = SY*VN(K) 

DIFF=VN(K)-YLS(K) 

100 bO-i>u+ ( Anr. ( D IM ))**2.0 

0K.MA=(SG/M«»(C.B)/(£Y/N) 

ME TcRN 
PND 


CUiiMJUTInC CmF AN( ImT I , IMF ) 

CCMMON /AAA/X(‘5l2) ,YN(£>12) .YLS(512).0 Y(‘jI2 ) ,N ■^M.S1oMA,<M'' * 

SYDE=0 ,0 

SYeoE=0.0 

KR -0 

2-10 continue " “ 

K R — K M ♦ I 
K =KM-KOa 1 

IF ( K.LT.2 , UK«KR,OT, I NT I ) GO TO 290 

5Y0E = SVDE^0. 5*( OY( K) ♦DY( K-l » ) * (X (K )-X( K-l ) ) " ■ 

SYEDE=SYEDE*0.5*(OV(K)aDY(K- 1) )«0«5*(X(K)4-X(K>n l*|X(K)-X(K-I ) ) 
GO TO 280 
290 CONTINUE 

KFr 0 -• -- • — 

3C0 continue 

kF =KF ♦ I 

k=km^kf 

IF (K.GT«N»0H.KF.GT«INTF» go to 310 “ - 

SYDE=SYDEAO. S*(DY(K» ♦OV(K-l) )*(X(K)-X(K-l) ) 

GYEUL=SYE0Ea0,5*(0Y<K1»DY(K-1) )*0«5*(X(K)«X(K-1) )*(X(K)-X(K.-1) ) 
GO TO 300 

3 10 C3NT INUE - ■ 

EM=SVEOE/SYOE 

HETUHN 

END 


SUl.wnUTINE SIMGl M.N.EPS. lER) 
C DIMENSION A(MM).R(NM) 

COMMON /S I M/ A I I 00 ) »R ( 1 0 ) 

K»' = N*M 

IF(M) 23.23.1 

1 IFH=0 
F I V=P. 

L J 1 L = 1 . M M 
TU - 4ES( A(L ) ) 

IF ( Ta-MI V ) 3.3,2 

2 PIV=TB 
I =L 

3 continue 

TOL =CPSAPI V 
LST =1 

OO 17 N = 1 , M 
IF (PIV) 23. 23. A 

4 IF ( ILR) 7.5.7 

5 IF (PIV-TCLl P.6.7 

6 ILR=K-l 

7 FIVI=1 ./A( I ) 

J=( I - I ) /M 

1 = I . j«m-K 

J= J + 1 -K 


• ^ ■ 
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I 


1 


I 


I 


CO e L =K .NM, y 
LL=L* I 

TH=PIV I •R(LL ) 

K( LL» =R( L » 

S fi(L)=TO 

IF (K-M) 9.1 'i. 16 

« LcNO=LSTtM-K 
IF( J) 12. 12. tC 

10 TI*J*M 

DO 11 L=L‘'.T.L£ND 
TQ= A( L ) 

LL=C>I I 
A(L ) = A (LL ) 

11 A(LL) =TB 

12 DO 13 L-LST,yy.M 
LL = l_ ♦ I 

TB = P I VI *A(LL ) 
A(LL)=A(L ) 

M A(L l = Tfl 
ACLST ) = J 
PIV=C «0 
LST=LST^1 
J = 0 

CO Ifc II=LST,LtND 
PIV I=-A( I I ) 

IST=1 I ♦M 
J = J 

CO 15 L=IST.MM.}« 

LL =L-J 

A(L ) = A(L ) ♦PIV I*A(LL) 
TB= ABS ( A(L ) ) 

. IF (Tb-PIV) 15.15.14 

14 PIV=TU 
I=L 

15 CONTINUr 

CO 16 L=K.NM.y 
LL=L* J 

16 R(LL ) =P (LL ) 4PI V I OR (L ) 

17 LST=LST+M 

18 IF <M-l) 23.22.19 

19 IST=MM4-M 
LST=M+1 

CO 21 1=2. M 

I I =LST- 

I ST=I ST-LST 

L=IST-M 

L = A ( L ) ♦ • 5 

CO *1 J=ll.^y.M 

Tb=R( J » 

LL = J 

DO 20 K=iST.yy,y 
LL=LL+ 1 

?0 TU = TB-A ( K ) (LL ) 

K = J + L 

R( J ) =P ( K ) 

21 R(K)=TB 

22 FETIJFN 

23 IER=-1 
RETURN 
END 
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B . RUN parameter 

DAI A CARD 

//I READ (NSET, KMON, KDAY, KYEAR, THETAI, THETAR, PHI) 

NSET (IIP) 

> 2: Run with calibration parameters 

(or-1) (Data Card # 2 ) 

■ 2: Run without re-calihrations 

KMON, KDAY, KYEAR (3110) 

MONTH /DAY /YEAR 

THETAI. THETAR. PHI (3F10.3) 

Angular parameters (8^.6^,4>) 


n 

(Req. if 
NSET > 2 
or ■ 1 


READ (El. M. NCVI. VI. NCV2. V2. NI, XLL. NF. XHL, 
INTF. KPG, KSTEP. NPW, NSTEP) 

El (F5.2) ; Incidemt-beam energy 

M(I5) : Order of the lest-square fitting program 

NCVI (13) 


Calibration of x-coordinate 



Operational limits 


Integration Limits 

INTF (15) 

KPG (15) ; Graphical Index 
> 1 : with graphical output 

< 1 : without graphical output 




INTI. 
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KSTEP (15) ; Step-size for graphical output 
NPW (15) : Output Index 

> 1 : with detail output 

< 1 : without detail output 

NSTE? (15) ; Step-size for detail output 
//3-#28 DATA SET: For thermal background gas 

#29-#5A DATA SET: For reflected beam 

ENDCARD: (Blank) 
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